Late time tails from momentarily stationary, compact initial data in Schwarzschild 

spacetimes 



Richard H. Price 

Department of Physics and Astronomy and Center for Gravitational Wave Astronomy, 
The University of Texas at Brownsville, Brownsville, Texas, 78520 

Lior M. Burko 

Department of Physics and Astronomy, Bates College, Lewiston, Maine O424O 

(Dated: February 5, 2008) 

Abstract 

An ^-pole perturbation in Schwarzschild spacetime generally falls off at late times t as t^'^^^'^ . It 
has recently been pointed out by Karkowski, Swierczyiiski and Malec, that for initial data that is 
of compact support, and is initially momentarily static, the late-time behavior is different, going 
as By considering the Laplace transforms of the fields, we show here why the momentarily 

stationary case is exceptional. We also explain, using a time-domain description, the special features 
of the time development in this exceptional case. 

PACS numbers: 04.70.Bw, 04.25. Nx, 04.30.Nk 



I. INTRODUCTION 

The perturbations in Schwarzschild spacetime radiate into the horizon and out to future null infinity (Scri-|-), so 
that at any fixed position in Schwarzschild coordinates, the perturbation falls off in time t. It has long been known 
that the fall-off is an inverse power-law in time of the form . For a multipole perturbation of multipole index H., 
the value of the power law index n can take several forms. In the case of a perturbation that, at large radius, has the 
asymptotic form of a static multipole, n = 2t + 2. This case is of particular astrophysical interest since it describes 
the fate of initial multipoles coupled to a star that has been stationary, but undergoes gravitational collapse to a black 
hole This result has been confirmed, for example, by Cunningham, Price, and Moncrief 0, and more recently by 
Baumgarte and Shapiro, in the context of the collapse of magnetized neutron stars Q. 

In the case of an initial moment that is asymptotically static, the initial field at large radius is the limiting factor 
in the rate at which the field falls off. The fall-off is faster if the initial data has compact support. The rule in this 
case IS n = This result was found numerically also for fully-nonlinear spherical collapse of a scalar field 

Recently, however, Karkowski, Swierczyiiski and Malec Q, hereafter KSM, presented numerical evidence that 
if the data are momentarily stationary as well as being of compact support, then n = 2i? + 4. Here we will explain, 
from two different points of view, why the momentarily static case is an exception. 

Perturbations of spherically symmetric black holes can be decomposed into multipoles, and each multipole moment 
^(x, i) satisfies an equation of the form 
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y(x)* = 0. (1) 



In the specific case of perturbations of a Schwarzschild spacetime of mass M, the variable t is the usual Schwarzschild 
time coordinate and x is the 'tortoise' coordinate, related to the Schwarzschild radial coordinate r by x = = 
r -|- 2Mln(r/2M — 1). Here we use units in which G = 1 = c; we choose M = 1 without loss of generality, so that 
t and X are dimensionless. If 4* represents odd-parity gravitational perturbations, then the potential V{x) is the 
Regge- Wheeler potential; if 5" represents even-parity gravitational perturbations, V{x) is the Zerilli potential. 
For even- or odd-parity electromagnetic perturbations, or for scalar perturbations, V{x) has a somewhat different 
form. It will be convenient here for us not to specify at the outset just what particular form V{x) takes. We will 
require only that V{x) falls off sharply as a; ^ —00 and that V{x) = £{£ + -|- 0[ln(a;)/a;]}, for large x. 

In the next section we work with the Laplace transform of '^{x,t), and relate the Laplace transform to an integral 
over the initial data. In principle, the form of the late-time tails of the perturbations can be extracted from the 
analytic details of Green functions in Laplace or Fourier space, as others have shown ■ Indeed, as shown by 

Leaver the late time tails will be a sum of a t~'^^~^ tail and a tail, the latter term arising from time- 

symmetric initial data. But we can avoid such complications. If one accepts that n = 2^ + 3 is the result for generic 
initial data of compact support, it turns out that an immediate consequence is that n must be 2£ + A if the initial 
data are momentarily stationary. 
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The late-time tails are usually thought of as a result of backscatter of radiation by the potential V{x) at large radius. 
Though the proof in Sec.^is definitive, it does not explain how, in the scattering picture, the momentarily stationary 
initial data are exceptional. In Sec. IIIII we provide a heuristic explanation by showing that for the momentarily 
stationary case, the initial data result in two outgoing pulses that, in a sense, cancel each other. 

The exceptional behavior of time-symmetric initial data can suggest the following paradox: Take the effective 
potential to be that of a Schwarzschild spacetime, but truncate it below a certain value of the Regge- Wheeler 'tortoise' 
coordinate , and take this truncation to be at a large negative value of r» . For the generation of tails, such a truncated 
potential is expected to be an excellent approximation to the full Schwarzschild potential, because the Schwarzschild 
potential drops off exponentially with for large and negative values of , and as is well known, it is only the form 
of the effective potential at large distances (large and positive values of r») that is important for the tails problem in 
Schwarzschild. Consider first initial data of an outgoing pulse of compact support to the "left" (more negative side) 
of the truncated potential, so that the initial pulse is fully located in the region of zero potential. One could expect 
the tail in this case to be given by t~'^^~^, since the initial outgoing pulse is generic time-asymmetric initial data. 
Consider next the same situation, but this time with an initially momentarily static pulse of twice the amplitude of 
the initially outgoing pulse we previously considered. The compact initial pulse is in a region of spacetime with zero 
potential, and therefore will immediately split into outgoing and incoming pulses. The latter is never heard from 
again; it travels to the left in a region of zero potential, and therefore never scatters. The pulse traveling to the right 
is identical to the situation we considered above. However, in this case, based on the prediction of KSM, the tails 
should be given by t~^^~'^. How can we explain this paradox, and what is the correct form of the tail in this situation? 
We conclude this paper by resolving this conflict of predictions. 



II. RELATION OF TAILS TO INITIAL DATA 



We now follow the approach used by several authors 0,0], and introduce the Laplace transform ■(/'(a;, s) of 'I'(x, t) 
through 



Jo 



(2) 



and the inverse 



C~^[i;{x,s)] = -^{x,t) 



2m 



e''*-^]{x, s) ds , 



(3) 



where F is a vertical contour in the right half of the complex s plane. With the relation C{d^/dt) = — 4'(a;,t 
0) -|- sipix, s), and its extension to second time derivatives, we write the Laplace transform of Eq. Q as 



^ - [s^ + Vix)] ^ = Six, s) - -5*0(2;) - «'o(a:) 



(4) 



Here *o and '^q are, respectively, the initial (t = 0) value of '^{x,t) and the initial value of d'^{x,t)/dt. 

To solve Eq. we again follow the approach of several authors we define homogeneous solutions yLix, s) 

and Unix, s) of Eq. Q that, respectively, represent waves moving inward through the horizon, and outward at spatial 
infinity: 



X — > — oo 



VL ~ e yB^ ^ e . (5) 

The Green function can be constructed in the usual way from and yji, and the solution to Eq. I^J is given by 

1 



^Pix,s) = 



Wis) 



ynix^s) yLix' ,s)Six' ,s)dx' + yLix,s) I yuix' , s)Six' , s) dx' 



(6) 



W{s) being the Wronskian determinant of the homogeneous solutions. We next use the form of S{x, s) in Eq. Q to 
write Eq. ^ as 



V'(x, s) = il)Qix, s) + ijjoix, s) , 



where 



tpoix, s) 
ipoix,s) 



Wis) 

-1 

Wis) 



yRix,s) / yL{x',s)'i/oix')dx' + yL{x,s) / y^ix' , s)-^oix') dx' 

J —00 J X 

j*x 1*00 

ynixjs) yL{x' ,s)i/oix') dx' + yL{x,s) ynix' , s)i/oix') dx' 



(7) 

(8) 
(9) 
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Now let us suppose that ^'o(a;) and *o(a^) are arbitrary (bounded) functions of compact support, and for every 
choice of these functions, except perhaps the choice vE'o(3^) = 0; let us suppose that the fields at any value of x fall off 



as t 



From this we conclude that for any bounded f{x) of compact support the expression 



F(x,s) 



-1 



W{s) 



yR{x,s) I yL{x\s)f{x')dx' + yL{x,s) I yR{x\ s)f{x') dx' 



(10) 



gives the Laplace transform of a function that falls off in time no slower than t^^^^^. Now note that sF{x, s) is the 
Laplace transform of the time derivative of this function, and that the time derivative will fall off as i"^^"**. We can 
therefore conclude that if the generic late time behavior is t~'^^~^ , then ipoi^, s) is the transform of a function that 
falls off as t~'^^~^ and ipoix, s) is the transform of a function that falls off as t^^^^"'. In the exceptional case that the 
initial data is momentarily stationary, ipoi^i *) vanishes, and the late time behavior is 



III. THE BACKSCATTER OF MOMENTARILY STATIONARY COMPACT INITIAL DATA 

To explain the n = 2£ + 4 tails we shall use the general heuristic framework developed in Refs. 0,^3: consider 
that there is a background problem, with a zero-order potential 

I _(£{£+ l)/x'^ x>xo , 
^0 x<xo ■ ^ ^ 

We will consider the remainder of the potential to be a perturbation, so that V{x) = Vq +e SV. The e is an accounting 
device so that we carry out a sort of perturbative analysis. 

The idea of this division of the potential into a background part and a perturbation is that the background is a 
pure centrifugal potential that cannot produce long-lived radiative tails. The tails, therefore, must be due to SV. 
The (scalar field) monopole case is somewhat awkward, since the centrifugal potential vanishes. What we really need 
though is some "edge" at some xq. Barack 'F] discusses the possibility of using a delta function for this purpose, but 
we need not be specific. Despite its awkward feature, we shall rely heavily on the monopole case. This is not solely 
because the slowly decaying £ — tails are the easiest to compute to long times. More important, the description 
of backscatter for the £ = case, lacks technical complications of higher order multipolcs. To focus on the essential 
ideas of backscatter, we shall therefore confine ourselves to £ — 0. The extension to higher £ is straightforward. 

We shall confine ourselves to descriptions to first order in e. This, sensu stricto, is not correct but we believe that 
the fundamental picture that comes out of that first-order analysis is correct. Strong evidence for this is the numerical 
accuracy (illustrated below in Fig. |2Jl of a prediction coming from this picture. 

As in Ref. 0, we shall introduce advanced time v and retarded time u by 

v = t + x u = t~x, (12) 

and we shall focus attention, not on tails at 1+ (that is, t — > oo at constant x), but rather at Scri-|- (that is, w — > oo 
at constant u). It can be shown P,[5l| that the tails at Scri-|- and at Z+ are tightly connected. If the former is m~™ 
then the latter is t~"^~^~^. For the monopole case, then, we need to show that for generic initial data of compact 
support, the tail at Scri-f falls off as it^^, while for momentarily stationary initial data of compact support the tail 
has the form u~^. 

To simplify some statements in our analysis, we will not deal with the monopole potential per se, but rather, shall 
take our potential to be strictly zero for x < xq, and to be 1/x^ for x > xq. We now write Eq. to first order in e 
as 

^iiluv = -\eSV^o- (13) 

Let us suppose that the zeroth order solution is an outgoing pulse 5* = F(u) of compact support. Following the steps 
in Ref. 5], and making the same approximations, we get 

= / F{u')du' dv'5V{x'). (14) 

^ J — OO J u-{-2xq 

When we use our special form of the potential 5V — this becomes 

1 PU, OO O /* OO 

= nu')du' ^-'t^t^"-^/ 

^J-oo Ju+2xn (V — U ) U J^ao 
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FIG. 1: The evolution of momentarily stationary initial data of compact support. The figure shows the zeroth order pulses 
Fi(w), G{v), and F2{u) described in the text. 

We have used the same approximations here as those of Ref. j^. In particular, we have assumed that the value of u 
to which this tail result is to be valid satisfies u S> u', where u' is any point in the support for F 14]. 

We now consider the case of initial data that are momentarily stationary. The zero-order solutions for such initial 
data will immediately "split" into an ingoing pulse and an outgoing pulse. These are labeled as Fi{u) and G{v) in 
Fig.n (Note that v has been conformally rescaled in the figure to bring Scri-|- to a finite location.) In the special case 
of momentarily stationary initial data the ingoing and outgoing zero-order pulses will be related by G{^) = Fi{—^). 
The "edge" at x — xq is a zero-order feature, so the ingoing G{v) will undergo partial reflection at xq and will generate 
a second outgoing pulse F2{u). 

For there to be no it~^ tail at Scri-|- (and hence no tail at X+) it must be the case that [-Pi(w) + -P2(u)] du — 
0. We have numerically checked a large number of examples, with different potentials, and different momentarily 
stationary initial data. In all cases we have found that the "cross section" (i.e. , the u integral) of the reflected pulse 
F2(u) is opposite in sign to Fi(u) and to numerical accuracy is equal in magnitude. Since G'(^) — fi(— this is 
equivalent to 



In practice, we integrated at t = const over the outgoing part of the field. To have a numerically zero integral, the 
reflected field F2 must have the opposite sign to the initial field Fi. In Fig. [21 we show the field at different values of 
the time as a function of the 'tortoise' coordinate. The figure shows the field soon after the time-symmetric initial 
data split into outgoing and incoming fields. As the outgoing field arrives at the peak of the effective potential the 
field scatters, and part is reflected toward the left with the opposite sign (and never heard from again), and a fleld of 
the opposite sign continues to move toward the right, following the main pulse. That is, the outgoing field is composed 
of the prompt field, and a broadened field of the opposite sign. It is the integral of the combined outgoing field that 
we calculate, and the result is shown in Fig. 21 In practice, we compute the integral only for positive values of r*, 
to capture only the contributions from the outgoing field. However, at and near r, = there is no clear separation 
between outgoing and incoming fields, and the field there is not strictly zero for finite values of time. Because of the 
contributions from the neighborhood of the peak of the effective potential, the integral does not vanish at finite values 
of time. However, the "area" between the field and the horizonal axis drops with time, and as t — > 00 the integral 
approaches zero, like 

This should not be misinterpreted as total reflection of the ingoing pulse G{v), in the sense of total reflection of 
energy in the waves. Such a statement about reflection refers to a quantity quadratic in the wave pulse; the "reflected" 
pulse F2(u) is generally quite different in shape from the ingoing pulse G(v), so Eq. (|16|l is very different from a claim 
of total energy reflection. 

The relationship in Eq. IjlGI) can be said to be the explanation for initially stationary initial data being a special 
case, and therefore of some importance. Though numerical veriflcation of this relationship is its ultimate justiflcation, 
it is interesting that there is a heuristic argument for Eq. (|16|) that helps us to understand it. Figure 01 shows a 
zeroth-order ingoing wave reflecting off the edge a,t x = xq. For this situation, let us integrate the relationship in 
Eq. ()13fl over the range ofu,v shown in Fig. 01 as the rhombus with vertices A,B,C,D. The right hand side is clearly 




(16) 
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FIG. 2: The field at four diff'erent value of the time. The initial data is that of a time-symmetric field centers at r, — — 20M, 
and the field is shown for t = 10 Af (dotted), t = 20M (dash-dotted), t = 40M (solid), and t = 60M (dashed curve). 




t/M 



FIG. 3: The "area" integral of the outgoing field vs. the time. The circles represent the numerically integrated area, and the 
solid reference line is TAi/t^. The errors in the numerical data are smaller than 1%. 



of order e. We make this explicit by writing 



5F*(o) dudv . 



(17) 



The integral on the left can immediately be evaluated: 

J y"*(i),™dwt^?^ = *(i)A + *(i)c-*(i)s-^'(i)D- (18) 

From causality we have 

*(i)c = = 0. (19) 

By taking point A at sufficiently large u, we can make = ^(i)^ arbitrarily small. Point B is is at /+. The scattering 
of the ingoing pulse to I"*" is zero, so = 0. 

We conclude that the left hand side of Eq. l(T7|l vanishes, and this means that the integral on the right of Eq. lfT7|) 
must vanish. We can break the right hand side integral into the contributions due to the zeroth order ingoing pulse 
G{v) and the zeroth order outgoing pulse F(u). For the ingoing pulse 



SVG{v) dudv ^2 
If we add the outgoing contribution, we find that 



dVG{v) dxdv^2 I dV{x) dx I G{v) dv . 



2 / 6Vix)dx 

' Xo 



G{v) dv + / F{u) du 



(20) 



(21) 
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FIG. 4: Reflection of an ingoing zeroth-order wave at the edge of the potential. 




FIG. 5: A spherical scalar field in Schwarzschild. The initial data are for an outgoing compact pulse centered at r«/M = —40. 
The solid line is the field, the dashed reference line is proportional to t~^, and the dotted reference line is proportional to t~'^ . 

must vanish, and hence we have given a heuristic explanation for Eq. (|16|l . 

We are now in a position to revisit the paradox we described in the Introduction. There is a fundamental, although 
subtle, difference between the case of the truncated potential and the true Schwarzschild potential: In the latter case 
there is a small overlap of the initial data and the rapidly decreasing effective potential. Because of this overlap, there 
is no exact equivalence between outgoing initial data and time-symmetric initial data of twice the amplitude. Because 
the overlap is small, we expect the field from an initially outgoing pulse to fall off as at intermediate times. 

However, the small overlap implies that the true late time behavior will be t^'^^"^ . This situation is demonstrated 
in Fig. |S1 which shows the results for an initially outgoing field that is located to the left of the peak of the effective 
potential, for a spherically symmetric scalar field. At intermediate times the field clearly falls off like i^"*, but the 
asymptotic fall off is , as expected. 
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